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Abstract 

In this paper, the conventional Holstein-Primakoff method is generalized with the 
help of the characteristic angle transformation [Lei Zhou and Ruibao Tao, J. Phys. A 
27 5599 (1994)] for the spin-one magnetic systems with single-ion anisotropies. We find 
that the weakness of the conventional method for such systems can be overcome by the 
new approach. Two models will be discussed to illuminate the main idea, which are 
the "easy-plane" and the "easy-axis" spin-one ferromagnet, respectively. Comparisons 
show that the current approach can give reasonable ground state properties for the 
magnetic system with "easy-plane" anisotropy though the conventional method never 
can, and can give a better representation than the conventional one for the magnetic 
system with "easy-axis" anisotropy though the latter is usually believed to be a good 
approximation in such case. Study of the easy-plane model shows that there is a phase 
transition induced by the external field, and the low-temperature specific heat may 
have a peak as the field reaches the critical value. 
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1 Introduction 



Magnetic systems with single-ion anisotropy D(S*) 2 have been attracting attentions 
for years since such kind of anisotropy was found to be very popular in many magnetic 
materials On theoretical side, the spin- wave excitation in such systems are 

not very easy to handle caused by the off-diagonal effect of the single-ion anisotropy 
especially when the spontaneous magnetized direction is not the same as the anisotropic 
direction. Many theoretical approaches have been developed to deal with such kind 
of systems. Usually, one first apply a rotating transformation of the spin vectors to 
determine the ground state, then perform a Holstein-Primakoff (H-P) transformation 
to study the low-lying spin wave excitations. However, this method was found to be a 
good approximation only when the anisotropy is the "easy-axis" case (i.e. D < 0). In 
the "easy-plane" case (i.e. D > 0), the method was much worse. To understand it, one 
can study an easy-plane Heisenberg model. If the conventional H-P method is used 
naively to discuss the ground state and the magnon excitations of such a system, an 
imaginary value of the excitation energy for u k = 0" mode will always be encountered, 
which implies the failure of this method. 

Such deficiency of the conventional method is caused by missing an important 
quantum effect. Actually, for the single- ion anisotropy (no matter "easy-axis" case or 
"easy-plane" case), off-diagonal term D sin 2 9(Sf ) 2 will always appear in the Hamil- 
tonian as well as the diagonal terms D cos 2 6(Sf ) 2 after introducing the spin vector 
rotation. Such off-diagonal term may have the tendency to mix the single-site spin- 
state \n) with \n + 2) and \n — 2) to form the proper eigenstates, and this spin-states 
mixing effect is completely a quantum one which is very important in "easy-plane" 
anisotropy case. Unfortunately, such a quantum effect has been neglected by the con- 
ventional H-P method. As the result, the conventional method failed for the magnetic 
systems with "easy-plane" anisotropy. 

On the other hand, many methods have been proposed for the easy-plane magnetic 
systems Pj-[TH]. The matching of the matrix elements (MME) method P|-[ffl] was one 



which can be used to consider the spin-states mixing effect perturbatively so that it 
can give a reasonable representation for an easy-plane ferromagnet when the single-ion 
anisotropy is small 0-0, and some numerical methods were developed for an easy- 
plane spin-one ferromagnet f|-[[l. Recently, another method - the characteristic angle 
(CA) method was proposed for the easy-plane spin-one ferromagnet which could be 
applied to describe such spin-states mixing effect by a variation parameter through a 



spin operator transformation [ID|. The magnetic properties had been investigated for 



such a system in zero field, and the results seemed to be closer to the numerical results 



than those of the MME method HTU 



The present work is focused on generalizing the conventional H-P method with the 
help of CA transformation for the spin-one magnetic systems with single-ion aniso- 
topies. Two particular models will be studied as the illustration of the CA approach, 
although the latter is certainly not limited to such models. The difficulties faced by 
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the conventional H-P method are overcome for such systems with the help of the new 
approach. 

This paper is organized as follows. In the next section, the easy-plane model is 
studied using the CA approach. Detailed comparisons of the CA approach with the 
conventional method is made in Section 3. Section 4 is devoted to an easy-axis model, 
and the conclusions are summarized in the last section. 

2 Easy-plane case 

The first model we will study is an easy-plane spin-one ferromagnet in an external 
magnetic field which is applied perpendicular to the "easy-plane" . The Hamiltonian of 
this system can be given as: 

H = - J £ S, . S, + D J2(S*) 2 - h £ St, (1) 

(hi) * * 

where the first term is the exchange interaction, and the second one is the single-ion 
anisotropy. The anisotropy parameter D is positive so that x-y plane is the so-called 
"easy-plane" and z axis is the "hard axis". An external magnetic field h is applied 
along the "hard axis" . 

Although the single-site part of Hamiltonian D(S*) 2 — hSf has already been a 
diagonalized form, it is still unreasonable to apply a H-P transformation naively to 
discuss the magnetic properties of such a system by assuming the ground state to be 
the ordinary ferromagnetic state. Actually if we do that, we will easily find that the 
magnon excitation energy of ll k = 0" mode will be always negative in the case of 
u h < D" . That is because the "starting point" based on which the spin deviations are 
discussed is wrong. 

One must be very careful in finding a reasonable "starting point" . Actually, in such 
a system, on the one hand, the spins are forced into the "easy-plane" by the single- 
ion anisotropy, on the other hand, they have the tendency to point along the "hard 
axis" caused by the external field. As the result, this two effects must compete with 
each other and a new direction z' axis would be optimized to describe the spontaneous 
magnetized direction. So, it's desirable to introduce a new coordinates system (x',y',z') 
in which the spin components are related to those in the original coordinates by the 
following transformation: 

S z t = cosOrSf -sm9 r Sf (2) 
Sf = cos r Sf + sin r Sf (3) 
S? = Sf (4) 
Applying the above transformation to Hamiltonian (1), we have 

H= - J XS'jiS'j + fl cos 2 9 r .J2(Sf) 2 + D sin 2 9 r . ]T(S 'f 'f -hcosO^Sf 

(i,j) i i i 
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- D sin 9 r cos 9 r J2( S f S f + S f S f) + h sin °r E S f ( 5 ) 

i i 

In a classical view, we can always determine 9 r based on variation method assuming 
all spins are aligned along the z' direction in the ground state. However, one should 
be careful in quantum case, especially in the current "easy-plane" anisotropy case. 
Actually, if we apply the H-P transformation naively to investigate the spin-waves 
excitation in such system, an imaginary value of the magnon excitation energy for 
u k = 0" mode will always exist. In fact, since 

D sin 2 9 r {Sff = ° sin 2 9 r (Sf Sf' + Sf' Sf) + | sin 2 9 r (Sf Sf' + SfSf'), (6) 

if the H-P transformation is applied naively to the Hamiltonian ([5]), one may find 
that the off-diagonal terms j sin 2 9 r (Sf Sf + Sf ) in the above equation have no 
contribution to the constant term of the transformed Hamiltonian. That means the 
spin-states mixing effect has already been neglected by the conventional H-P method. 
Unfortunately, such effect is very important and must be considered in such case. The 
characteristic angle (CA) transformation [0] was developed to describe the spin-state 
mixing effect in spin-one case by introducing another variation parameter 6 C : 

Sf = cos 6 c Sj~ + sin 6 c Sfexp(inSj), (7) 
Sf = cos6 c Sf + sin 6 c exp(— in Sj)Sj~ , (8) 
5/ = (l/2)[Sf,Sf').. (9) 

The spin operators are transformed to a new set of quasi-spin operators (Sp, S?) which 
had been proved to obey all spin-one operator's commutation rules ||10|| . After the 
CA transformation, we can apply a H-P transformation to transform the quasi-spin 
operator to Bose one, 

Sf — ► l-afa h (10) 
Sf — > V2Jl - (atai/2) ai , (11) 



Sf — v^a+Vl - (a^fli/2). (12) 
Then, the Hamiltonian will have the following form: 

H = U + H 1 + H 2 + ■ ■ ■ (13) 

where 

U = N{-JZ cos 2 29 c + D - ^sin 2 6 r (l + sin29 c ) - hcos6 r .cos29 c ] (14) 
y/2 

Hi = — — \D sin 9 r cos 9 r (cos 9 C + sin 9 C ) — h sin 9 r (cos 9 C — sin 9 C ) ] 

i 

(af + af (15) 
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and H2 can be written in momentum k space as follows 

H 2 = J2 A k a t a k + J2 B k( a k a -k + a ka-k) (16) 

k k 

A k = 2JZ(cos 2 9 c - 7 fc ) + y sin 2 9 r (l + sin 29 c ) + h cos r cos26> c 

-Dcos 2 9 r , (17) 

= [— JZ sin 40 c + — sin 2 9 r cos 29 c — h cos r sin 29 c ] 

+ JZsin29 c j k . (18) 

Based on the variation method we understand that the two parameters 9 r , 9 C should 
be determined by minimizing the ground state energy. As a first order approximation, 
we may obtain: 

■Ti^rUo(9 r ,9 c ) = -£>sin0 r cos0 r (l + sin20 c ) + /i sin r cos 20 c = (19) 
N d9 r 

T7~nr U o(8r,Q c ) = 4 JZ sin 20 c cos 29 c - D sin 2 9 r cos 29 c + 2h cos 9 r sin 20 c 

= (20) 

Eq. (^) is just the same as the condition if 1 — and Eq. (|20|) can cancel most 
of the off-diagonal terms which are in the square bracket in the expression of B k . If 
we substitute the solution of the above non-linear equations into to the Hamiltonian 
(|13"D, then diagonalize the harmonic part of Hamiltonian H2 by a usual Bogolyubove 
transformation, the total Hamiltonian will be: 

H = U^ + J2E k aia k + --- (21) 

k 

where 



U'o = Uo- 1 -^ A k + \Y,^l-ABl (22) 

k k 



E k = y/Al-4Bl (23) 

The ground state in such a method can be defined by 

a*|0) = (24) 

Then the induced magnetization M(h) is derived in the harmonic approximation as 
follows: 



M(h) = ^E(oi^|o)^^E c ^^(o|sf|o) 

i 
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— 5>os0 r {cos20 c (O|£f|O) + sin# c cosfl c (0|(,S+) 2 + (S^) 2 \0)} 



— ^cosfl r {cos2# c (0|l - afai\0) + V2 sin6* c cos6 l c (0|a^ 2 + af |0) } 

i 

3 1 ^ cos # r cos20 c A fe + 2y / 2cos # r sin 26> C I?A; 

- cos 9 r cos 20 c > = . (25) 

2 2iv * v^r^! 



Thus, put the solution of Eqs. ©-(f^) into Eqs. (g|)-(gD, such physical 

properties as the ground state energy, the magnon dispersion relation and the induced 
magnetization can be obtained. However, since it is very difficult to solve the non- 
linear equations analytically, numerical calculations are carried out. The system with 
anisotropy parameter D/4JZ = 0.6 has been studied as an example. 

9 r and 9 C as the functions of the external field have been drawn together in figure 1 
from which one can find that they are both the decreasing function of the external field. 
It is understood that 9 r is used to describe the spontaneous magnetized direction and 
9 C the spin-states mixing effect, in zero applied field case, the spontaneous magnetized 
direction will be the x axis (9 r = 90°) and the spin-states mixing effect should be the 
strongest since the off-diagonal term D sin 2 9 r (Sf ) is the strongest, and the value of 9 C 



is consistent with Ref. [ID where h = case has already been discussed. While the 
external magnetic field is strengthened, the spins will point along a direction which is 
closer to the z axis due to the interaction with the external field so that 9 r will decrease, 
at the same time, the spin-states mixing effect is also weakened since the off-diagonal 
interactions in the total Hamiltonian will turn smaller along with the decrement of 
9 r . However, when h reaches a critical value h c = D, the external magnetic field is so 
strong that the spins will not be rotated any longer, and the off-diagonal term comes 
to zero either, as the result, 9 r and 9 C will vanish simultaneously. 



3 Comparisons and Discussions 

In this section, we will compare the CA method with the conventional H-P method in 
details and discuss the magnetic properties of the above mentioned system. 



First, one may find that more quantum effects have been comprised in the constant 
term of the Hamiltonian by the new approach. 

Introducing the H-P transformation naively to Hamiltonian (5), the constant term 
can be found as: 

U^ p = N(-JZ + D cos 2 9 r -hcos9 r )+ N— sin 2 9 r 

= Uq + N^- sin 2 9 r . (26) 
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where Uq is the ground state energy obtained by a classical rotating transformation. 

After applying the CA transformation, the ground state energy is Eq. (|14]) which 
can be rewritten as: 

U = U^ P + U x (27) 

where 

C/i = N(JZsm 2 29 c - ^ sin 2 fl r sin 2# c + 2/i cosfl r sin 2 #). (28) 
Ui is an additional term introduced by the CA transformation which will vanish as 

e c = o. 

The conventional H-P method is a semi-classical one. The only quantum effect in 
fl26|) is the term iVy sin 2 9 r which comes from the contribution of sin 2 9 r (Sf + 
S~ Sj~ ) in Eq. (|G|), and other terms which have been collected in Uq of Eq. ( P5| ) can be 
easily recovered by a classical method. However, after applying the CA transformation, 
it can be clearly found that there is an additional contribution U\ in the expression 
of Uq which describes the single-site spin-states mixing effect through the variation 
parameter 9 C . Such an effect is completely a quantum one which has no classical 
counterpart, and it is expressed as the competition of the exchange term (JZ) and the 
external term (h) with the single-ion anisotropy term (D). 

So, more quantum effects have been considered by the CA approach than the con- 
ventional H-P method even in the constant term of the Hamiltonian. Furthermore, 
considering such quantum effect will lead to a lower ground state energy since the 
ground state in CA method is selected to be the minimum point of function Uq al- 
though that in the conventional H-P method is not so. 



Now, one can compare the CA method with the conventional H-P method for the 
elementary excitation of such system. Put the solution of 6 r and 9 C into Eq. (|23|) , 
the magnon excitation gap can be calculated with respect to the external magnetic 
field and the result has been shown in figure 2, where one can find that the magnon 
excitation gap obtained by the CA method will always be positive or zero. 

However, based on the conventional H-P method, one should obtain the variation 
parameter 9 r by minimizing (|2q) which yields: 

4r U o HP = ~ D sin °r cos 6 r + h sin 9 r = (29) 



then substitute the solution of the above equation into Eqs. ( |T7HT8| ), (p3|), one can 



easily find that the magnon excitation gap in H-P method will be 



AffP h 2 -D 2 , . 

= V 2D ' (30) 
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Of cause, the excitation gap will never be real when h < D. That is to say the 
conventional H-P method can not be applied naively to study the magnetic systems 
with "easy-plane" anisotropy. However, the CA method has overcome this difficulty as 
shown in figure 2. 



The induced magnetization as the function of the external magnetic field has been 
drawn in figure 3. From figures 1-3, one may find that the point h c = D is very 
strange and there seems to be a phase transition in such a point. As the external field 
is strengthened across h c , the system transits to a phase in which the spins are not 
rotated any longer. Actually, this phase transition can be clearly shown by calculating 
the low-temperature specific heat C v . 

Suppose the system is at low temperature, and only the low energy excitation is 
considered, then the inner energy will be: 



So the specific heat can be obtained: 



dE(T) . 1 ^ (j^f ) 2 exp( 



E k 



The specific heat of the system as the function of the external magnetic field at the 
temperature KbT/JZ = 0.1 has been shown in figure 4, in which a peak can be 
apparently found in the critical point h c = D. The physics can be understood as 
follows: in the vicinity of the phase transition point h c , the magnons will be excited 
without a gap (fig. 2) so that the thermal fluctuations are strong. 



4 Easy-axis case 

Now we will study another model which is an "easy-axis" spin-one ferromagnet in 
an external magnetic field whose direction is perpendicular to the "easy-axis". The 
Hamiltonian of such system can be given as: 

H = - J £ Si • S, - D ) 2 - h £ St, (33) 

i i 

where the single-ion anisotropy makes x axis an easy-axis, and an external magnetic 
field is applied along z direction. 

At a first glance, this Hamiltonian has a same classical picture as the last model - 
the anisotropic interaction and the interaction with external field have to compete with 
each other and will be balanced at some angle 9 r . So a rotating transformation of the 
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spin vectors is helpful. In fact, many authors have used this method to discuss various 
kinds of magnetic systems with single-ion anisotropy and believed this approximation 
will work well when the single-ion anisotropy is "easy-axis" case. After the rotating 
transformation (§)-(§), the Hamiltonian will be: 

H= J J2 S' \ «S' 'j - Dsm 2 9 r Y,{S* 'f-D cos 2 r ^(Sf ) 2 - /i cos OrY.s? 

i i i 

- D sin 9 r cos 9 r ^2(SfSf + SfSf ) + h sin 9 r £ Sf' (34) 

i i 

there are still off-diagonal terms (— D cos 2 9 r J2i(Sf) 2 ), in the Hamiltonian, and this 
off-diagonal interactions may be important in some cases. So it is helpful to apply the 
CA transformation to get a more reasonable representation. 

Actually, after almost the same procedure as that for the easy-plane model, the 
Hamiltonian can be transformed to: 

H = U + H 1 + H 2 + ■ ■ ■ (35) 

where 

U = N[- JZ cos 2 29 C -D + — cos 2 9 r (l + sm29 c )-hcos9 r cos29 c ] (36) 

2 

Hi = y~l[~ D sinfl r cos^ r (cos6 l c + sinfl e ) + fcsinfl r (cosfl c — sin^ e )] 

{at + <0 (37) 

and 

H 2 = ^2A k a^a k + ^2B k (a^at k + a k a- k ), (38) 

k k 

A k = 2JZ{cos 2 29 c — 7^) — — cos 2 9 r (l + sin 29 c ) + h cos 9 r cos 29 c 

+ Dsm 2 9 r , (39) 

\p2 D 
B k = —^-[—JZsm4:9 c — — cos 2 9 r cos 29 c — h cos 9 r sin 29 c ] 

+ JZ sin29 c j k . (40) 
The two variational parameters 9 r , 9 C satisfy: 

^-7^U (9 r ,9 c ) = -Dsm9 r cos9 r (l + sm29 c ) + hsm9 r cos29 c = (41) 
N d9 r 

^—U (9 r ,9 c ) = 4 JZ sin 29 c cos 29 c + D cos 2 9 r cos 29 c + 2h cos 9 r sin 29 c 
N d9 c 

= (42) 
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where Eq. (fH|) will cancel the Hi part of the Hamiltonian. 

Physical properties such as magnon excitation, the induced magnetization have 
the same forms as which in last model (i.e. Eqs. ( p2[]23]) ,(p5|)) except the concrete 
expression of the functions Ak and Bk- 

Very similar to the "easy-plane" case, the constant term in the Hamiltonian can 
further be divided into two two terms: 



U Q = + U x (43) 



where 



U* p = N(-JZ - Dsm 2 6 r -hcos6 r )-NjCos 2 6 r 

= - iVy cos 2 9 r (44) 
U x = N(JZ sin 2 29 c + — cos 2 6 r sin26 c + 2hcos6 r sin 2 6 C ). (45) 

Uq F is the contribution of the conventional H-P method, and U\ is an additional 
term which describes the quantum effect of spin-states mixing in single-site. All the 
discussions are similar to the first model: more quantum effects have been comprised in 
the constant term of the Hamiltonian by the CA method, and considering such effect 
in CA method will lead to a lower ground state energy than that in the conventional 
H-P method. Actually, as shown in figure 5 where Uq and U^ p are drawn together 
with respect to the external magnetic field for an easy- axis spin-one ferromagnet, Uq 
is always found to be lower than Uq P . 

Now it is interesting to compare the elementary excitations calculated by the CA 
method with those by the conventional H-P method. In the latter case, 9c = and 9 r 
is obtained by minimizing Uq P (|44|) which yields: 

TflT U o P = ~ D sin 9r cos dr + h sin 9r = ( 46 ) 

N UUr 



So, substitute the solution of the above equation back into Eqs. ([^-|4T]) then into 
Eq. (|23|), the elementary excitation in the conventional H-P method can be calculated 
readily. 

The magnon excitation gaps have been calculated by both methods with respect 
to the external magnetic field, and the results are presented in figure 6 where the solid 
line is by CA method and the dotted line is by the conventional H-P method. From the 
figure one may find that: when the external field is close to the anisotropy parameter 
D, there is a small region where the magnon excitation gap calculated by the H-P 
method will be imaginary, which indicates that this approximation is poor in such 
area. However, the solid line in the figure tells us that CA method has overcome this 
difficulty and the magnon excitation gap will always be real and positive in CA method. 
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Actually, as shown in figure 7 where the values of the two variation parameters 9 C 9 r are 
drawn together with respect to the external field, when h is close to D, 6 r comes to zero 
and 6 C becomes somewhat larger indicating that the spin-state mixing effect caused by 
the off-diagonal interaction may be very strong. So, we must consider such effect with 
the help of CA transformation in such case, otherwise, the "starting point" may be 
unreasonable and will lead to an imaginary minimum excitation energy. Outside this 
region, the off-diagonal terms are not so strong comparing to the diagonal parts, as 
the result, the spin-states mixing effect is not very drastic and the conventional H-P 
method might be a reasonable approximation as many authors believed. However, the 
CA transformation may always be helpful to get a more reasonable representation for 
such a system. 

5 Conclusions 

To summarize, in this paper, the conventional method has been generalized with the 
help of the characteristic angle transformation for the spin-one magnetic systems. The 
difficulties faced by the conventional H-P method for magnetic systems with single-ion 
anisotropy has been overcome by the new approach. Two models have been discussed 
to illuminate the main ideas, of which one is an easy-plane spin-one ferromagnet in 
an external field applied perpendicular to the "easy-plane", and the other is an easy- 
axis spin-one ferromagnet in an external field applied perpendicular to the "easy-axis" . 
Comparisons between the new approach and the old one show: more quantum effects 
have been considered by the CA method, as the result, CA method can examine the 
ground state properties of the "easy-plane" spin-one ferromagnet although the old 
method never can, and CA method can give an improved representation for the "easy- 
axis" spin-one ferromagnet although the conventional H-P method is usually believed 
to be a good approximation in such case. Also, study of the easy-plane model shows 
that a phase transition may take place induced by the applied field, and the low- 
temperature specific heat is found to have a peak when the external field reaches the 
critical value. 
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Captions: 



Figure 1: 6 r , 6 C as the functions of the external field for the easy-plane spin-one ferro- 
magnet with D/AJZ = 0.6. 

Figure 2: Magnon excitation gap as the function of the external field for the easy-plane 
spin-one ferromagnet with D/AJZ = 0.6. 

Figure 3: Induced magnetization as the function of the external field for the easy-plane 
spin-one ferromagnet with D/AJZ = 0.6. 

Figure 4: Specific heat as the function of the external field at the temperature K B T/ JZ = 
0.1 for the easy-plane spin-one ferromagnet with D/AJZ = 0.6. 

Figure 5: Ground state energy with respect to the external magnetic field in the case 
of using the conventional H-P method U^ p and using the CA method for an easy-axis 
spin-one ferromagnet with anisotropic parameter D/AJZ = 0.3. 

Figure 6: Magnon excitation gaps of the easy-axis spin-one ferromagnet with D/AJZ = 
0.3 as the functions of the external magnetic field for the case using the CA method 
(solid line) and using the conventional H-P method (dotted line). 

Figure 7: 9 r , 9 C as the functions of the external field for the easy- axis spin-one ferro- 
magnet with D/AJZ = 0.3. 
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